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Synchrotron-based x-ray reflectivity is increasingly employed as an in situ probe of surface mor- 
phology during thin film growth, but complete interpretation of the results requires modeling the 
growth process. Many models have been developed and employed for this purpose, yet no detailed, 
comparative studies of their scope and accuracy exists in the literature. Using experimental data 
obtained from hyperthermal deposition of pentane and diindenoperylene (DIP) on SiC>2, we com- 
pare and contrast three such models, both with each other and with detailed characterization of the 
surface morphology using ex-situ atomic force microscopy (AFM). These two systems each exhibit 
particular phenomena of broader interest: pentacene/SiCb exhibits a rapid transition from rough to 
smooth growth. DIP/Si02, under the conditions employed here, exhibits growth rate acceleration 
due to a different sticking probability between the substrate and film. In general, independent of 
which model is used, we find good agreement between the surface morphology obtained from fits 
to the in situx-ray data with the actual morphology at early times. This agreement deteriorates 
at later time, once the root-mean squared (rms) film roughness exceeds about 1 ML. A second 
observation is that, because layer coverages are under-determined by the evolution of a single point 
on the reflectivity curve, we find that the best fits to reflectivity data — corresponding to the low- 
est values of xl — do not necessarily yield the best agreement between simulated and measured 
surface morphologies. Instead, it appears critical that the model reproduce all local extrema in the 
data. In addition to showing that layer morphologies can be extracted from a minimal set of data, 
the methodology established here provides a basis for improving models of multilayer growth by 
comparison to real systems. 

PACS numbers: 68.43.Mn, 68.55.-a,81.05.Fb,81.15.Aa,83.85.Hf 



I. INTRODUCTION 

In situ, surface-sensitive scattering techniques, such as 
reflection high energy electron diffraction (RHEED) and 
x-ray scattering at the so-called "anti-Bragg" positiorfS^, 
yield direct information about surface morphology during 
growth, and have been applied to virtually all methods 
of thin film growth-such as electrodepositiorP, molec- 
ular b eam epitaxy (MBE), chemical vapor deposition 
(CVD) 5116 , and pulsed laser deposition (PLD^El Tne ir 
principal advantages, compared to scan-probe micro- 
scopies such as atomic force microscopy (AFM) , are their 
time resolution and the ease with which they can be in- 
corporated into growth systems as qualitative, in situ 
probes. For example in homoepitaxy, growth is classi- 
fied as "step-flow" (SF), "layer-by-layer" (LBL) or "three- 
dimensional" based, respectively, on whether the specu- 
lar reflectivity remains constant, oscillates, or decreases 
monotonically with increasing film thickness. An addi- 
tional advantage of x-ray scattering is the accuracy with 
which precise calculations of scattered intensity can be 
performed. Increasingly, this has been exploitecP^H to 
extract quantitative information about the surface mor- 
phology, namely the layer coverages as a function of time 
n (t), during growth. This analysis, however, requires a 
specific model of the morphological evolution of the film. 

Figure [l] shows x-ray reflectivity (XRR) and AFM data 
obtained from four pentacene films, grown in immediate 



succession on a single SiC>2 substrate under nominally 
identical conditions (substrate temperature 48° C), and 
serves to illustrate the main subject of this work: how 
accurately can the height probability distribution, which 
may be obtained directly from the AFM data in Fig.JlB- 
E, be extracted from fits to the x-ray data in Fig. 
In addition, how does the choice of model used for such 
analysis affect the simulated morphology? 

The models described here are necessary, in part, be- 
cause x-ray intensity measurements such as those in Fig- 
ure are performed at only one point in reciprocal 
space. As has recently been demonstrated-^ 7 * 18 * 22 * 23 *, 
measuring multiple points along the reflectivity curve can 
reduce or potentially obviate the requirement for such 
models. However, such measurements are not always 
possible or compatible with a particular system or experi- 
ment. Further, with notable exceptions^, such enhanced 
measurements are often used not to replace the use of 
growth models, but to limit uncertainties associated with 
comparing measurements to calculations. Finally, al- 
though this manuscript deals with the limited problem 
of determining the surface morphology of a growing film, 
the models described here also concern the broader chal- 
lenge of understanding, in detail, how microscopic, ki- 
netic parameters of a system determine its morphology. 
Clearly, comparing the simulated morphologies of these 
models to real systems is a critical, first step in validating 
them for deeper, more general use. 

One motivation for this work is to evaluate a par- 
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FIG. 1: (Color online) (A) XRR data obtained near the 
pentacene anti-Bragg position during pentacene deposition 
on Si02- Data are shown for four separate films grown to 
thickness of (A) 0.34 ML, (□) 1.4 ML, (x) 2.5 ML , and (+) 
4.4 ML, performed in immediate succession on the same sub- 
strate under the same conditions. (B-E) AFM images of the 
four films represented in (A) in order of increasing thickness. 
All images are 20 pimx20/im. 



For organic molecules containing only light elements, this 
method is not possible. Growth rate acceleration may 
also b e determined using ex situ microscopy (see section 
IV B| below) , but this approach requires several samples 
per growth condition, and can be severely complicated 
by the presence of de-wetting, a common phenomena in 
organic thin films^. Accurate determination of growth 
rates and growth rate acceleration from time-resolved 
XRR data relies directly on quantitative analysis as de- 
scribed above, and thus on the accuracy of the growth 
simulation. In prior studies, we validated our model by 
comparing both the total film thickness and rms rough- 
nesses of thin films measured by AFM to those predicted 
by best-fit simulations. Here, we provide a complete de- 
scription of this model, and also address the specific prob- 
lem of measuring growth rate acceleration. 

We proceed, in section [ill A| by reviewing XRR inten- 
sity calculations relevant for thin film growth. To eluci- 
date how different features of the film and substrate in- 
fluence the character of oscillations, we explicitly derive 
x-ray scattering parameters in the approximation where 
each layer consists of uniform density slabs. Additional 
insight is gained by considering the case of continuous, or 
step-flow growth, although it is not directly relevant for 
the particular systems described here. In sect ion [Tl I B| wc 
provide detailed descriptions of three distinct determin- 
istic, rate-equation models of thin film growth. Finally, 
in section [TV] we compare these models to two particular 
systems. The first of these, pentacene/Si02 (represented 
in Fig. [I]), has well-known structural and morphological 
evolu tion and, under these growth conditions (see section 
IV A[ ), does not exhibit growth rate acceleration. Four 
films of varying film thickness are grown under identical 
conditions. Height distributions of these films, measured 
with ex situ AFM, are then compared in detail with those 
deduced using fits to reflectivity data at the anti-Bragg 
position. The second system described in this paper, 
diindenoperylene (DIP)/Si02, exhibits growth rate ac- 
celeration, which we model as a difference in molecular 
sticking probability between the substrate and deposited 
film. We find that two of the models, including that in 
Ref. HH1 accurately describe the degree of acceleration. 



ticular model (the "modified Cohen" model, see sec- 
tion IIIB below) we recently developed to measure the 



time-dependent growth rate of hyperthermally-deposited 
organic thin films^. This time-dependence arises from 
differences between the sticking probabilities of de- 
posited molecules on the bare substrate and the grow- 
ing film. Because this sticking probability depends 
strongly on the substrate, a quartz crystal monitor is 
an inappropriate measure of growth rate. This phe- 
nomena has been observed in other systems: in metal- 
organic based GaN growth on sap phire, both growth rate 
accelerations! and deceleration^! have been observed us- 
ing in sziitmeasurements of Ga Ka x-ray fluorescence. 



II. EXPERIMENTAL 

Experiments were conducted at the G3 station at the 
Cornell Higher Energy Synchrotron Source (CHESS) in 
a custom-designed, ultrahigh vacuum chamber with a 
base pressure of 5xl0~ 9 mbar described previousl y 19 * 25 !. 
Growth was performed on 300 nm-thick Si02 thermal ox- 
ide films on Si (001) wafers, using a hyperthermal molecu- 
lar source. The deposition energy, measured using time- 
of-flight mass spectrometry, was approximately 2.5 eV for 
pentacene and 4.2 eV for DIP. Synthetic W/B4C multi- 
layer monochromators set the beam energy to 9.75 keV 
with AE/E 1.5%. The unattenuated flux at the 
sample position was approximately 5xl0 13 photons-s -1 - 
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-mm -2 . In order to avoid radiation damage, the beam 
was attenuated by a factor of 10 or more using multiple 
layers of aluminum foil at the upstream end of the hutch. 
An avalanche photodiode detector (Oxford Danfysik, Ox- 
ford, UK) scintillator counter was used for measuring the 
scattered x-ray intensity. AFM was conducted ex situ in 
tapping mode using a Digital Instruments 3100 Dimen- 
sion microscope (Santa Barbara, CA). 

III. THEORY 

A. X-ray Scattering 

The computational simplicity of x-ray scattering re- 
sults from the fact that, far from a Bragg peak or from 
the critical angle for total external reflection, the Born 
approximation (also referred to as the single scattering 
or kinematic approximation) may be usecpSJ. In this ap- 
proximation, the scattered intensity I(q) is proportional 
to |j4(g)| 2 , where the scattering amplitude A(q) is the 
Fourier transform of the electron density: 

Mo)= / p(r)exp(-iq-r). (1) 
Jv 

For the particular case of in situdata collected during 
thin film growth at a position q = q z z in reciprocal space, 
and with the assumption that the substrate does not vary 
during deposition, this may be simplified to l 14 l 16 l 23 l ; 

A(q z ,t) = A suh e 1 ^ + A mm £ 6 n {t)e-^< n - l \ (2) 

n=l 

Here, 6 n is the fractional coverage of layer n, and c 
is the lattice parameter of the film normal to the sur- 
face. Typically, experiments are conducted at the anti- 
Bragg position, corresponding to q z c = ir. Since Afii m 
and A su b are both defined as pure real, there are only 
two unknown scattering parameters, </> su b and the ra- 
tio ^4fiim/^4sub- These parameters may be calculated 
directly if the atomic structure-including the detailed 
structure of the interface-is known. For crystalline sub- 
strates, A au \ye l ^ sub incorporates a semi-in finite sum over 
the buried substrate lattice. Krause et a/! 27 1 28 1 explicitly 
calculated A(q z ,t) for the case of the organic molecule 
growth on Ag(lll). Alternatively, they can be treated 
as free parameters. 

In appendix [A] we derive a general expression for ^4 su b 
and (ftsub under the assumption that the sample geom- 
etry resembles that of Fig. [2j Specifically, we assume 
that the substrate is composed of a thin amorphous layer 
with density p± and thickness T\ — ct on an amorphous, 
semi- infinite layer with density po- This model is partic- 
ularly suitable for amorphous substrates such as Si02, 
but also valid for Si substrates with a thin Si02 na- 
tive oxide^, and thick, th ermal Si02 layers covered by 

a se lf-ass embled monolaye r^ 20 * 30 ! or an interfacial water 
layepM] 
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FIG. 2: Real-space schematic representation of organic thin 
film growth. Parameters c, r, A, n , and a n are described 
in the text. ki,kf, and q correspond to the incident and final 
wave vectors, and momentum transfer, respectively, in spec- 
ular geometry. 

To help visualize the scattering intensity arising from 
Eq[2]near the anti-Bragg position for organic thin films, 
it is conceptually useful to further approximate the film 
as composed of uniform density slabs of height c and den- 
sity P2- This approximation explicitly eliminates Bragg 
Peaks, but is nevertheless reasonably accurate near the 
anti-Bragg position whenever this geometry probes a 
length scale, I — 2Ti/q z = 2c, that is large compared to 
interatomic distances. For example, for pentacene and 
DIP, c sal. 54 nm and 1.66 nnPH, respectively. In this 
approximation, we have the result 

A suh e^ = -^-ie^ 1+2A \ Pl + (po - Pi)e i27rir ). (3) 
ZirLt 

0sub is simply the complex phase of A su i, above, while 
^4fiim/^4sub is 

Ajum _ 2/92 sin irL 

A suh ~ \ Pl + (po - pi)e*^| ' U 

In Eqs. [3] and |4j L = q z /(2n/c), and A corresponds to 
a small difference, as a fraction of unit cell height c, in 
height of the center of the first deposited layer from c/2 
above the substrate surface. Despite the fact that these 
equations include more than the two original parameters 
in Eq. [2j this representation generally reduces parameter 
space, since these parameters are at least approximately 
known, and can thus be constrained. 

Recently, Kowarik et al^l described the varying ap- 
pearance of anti-Bragg intensity oscillations caused by 
variations in the two parameters </> su b and Afii m /A su b. 
Equations [3] and [3] connect these parameters with details 
of the system. To help interpret and illustrate Eq. [2j the 
evolution of both the scattering amplitude and intensity 
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under four distinct scattering conditions are shown in 
Fig. [3] In all eight plots, the solid curve corresponds 
to the same evolution of coverages 9 n (t): namely, those 
obtained from a fit of the "modified Cohen" model (see 
section III B I to the data in Fig. [T]A. The dashed curves 
correspond to ideal LBL growth, while the dotted curves 
correspond to ideal SF growth. Although SF growth is 
not relevant for the examples discussed below, we nev- 
ertheless plot it to explicitly demonstrate the relation- 
ship between "roughness" and "Kiessig" oscillations. The 
scattering amplitude during SF growth is derived in ap- 
pendix [B] 

In each of the left-hand plots in Fig. [3J A su be l ^ sub is the 
vector originating from the origin, while A&i m originates 
at the tip of ^4 su be^ BUb and extends along the positive 
real axis. Figure [3]^ corresponds to the precise parame- 
ters obtained from the fit, as described below and shown 
in TablelU Because the best-fit value of L is somewhat less 
than 0.5, the scattering amplitude of ideal LBL growth 
traverses a non-convex regular polygon. The subsequent 
pairs of plots simulate other cases of interest. For exam- 
ple, Fig. [3j3 corresponds to the familiar case of homoepi- 
taxy, which exhibits two one oscillation per layer in LBL 
growth, and constant intensity during SF growth. Fig- 
ure [3p corresponds to growth on a 1 nm water layer 
on SiC>2. In this case, the Si02 substrate and water 
layer nearly cancel (po/pi ~ 2.2 =^ pi « po — p%, and 
exp(z27riT) « —1), so that the ratio Afn m /A su b becomes 
large (see Eq. |4| , with the result that the peak intensity 
near completion of odd-numbered layers is large com- 
pared to the starting intensity. Finally, Figj3]D simulates 
growth on SiC-2 as observed at the quarter-Bragg position 
(£ = 0.25). 

A clear feature of the solid curves in the left-hand plots 
in Fig. [3] is that the scattering amplitude of a rough film 
approaches the center of the circle corresponding to ideal 
SF growth. For L = 0.5, it has been shown previously^ 
that the scattered intensity of a rough film coincides pre- 
cisely with that of a film with precisely 0.5 ML coverage, 
i.e. #1 = 0.5, 9 n >i = 0. Figure [3] shows that this coinci- 
dence does not occur for L ^ 0.5. A more general inter- 
pretation of this limit is to view the scattering as arising 
from the different interfaces in the systerrP^ When the 
film/vacuum interface becomes rough, its effective scat- 
tering strength drops. Explicitly, evaluating Eq. |A4| for 
a rough film (e.g. taking 9 n — e~ Sn ) leads to 
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FIG. 3: (Color online) (left) Complex scattering amplitude 
and (right) intensity under identical growth conditions but 
different scattering conditions. In all plots, the dashed and 
dotted curves (color online) correspond to ideal LBL and 
SF growth, respectively, while the solid curve corresponds to 
growth according to the best-fit curve of the data in Fig.[TJ\. 
The four pairs of plots simulate (A) pentacene/Si02 identi- 
cal to that found from the fit to the modified Cohen model 
shown in Fig. [4] (B) homoepitaxy at the anti-Bragg po- 
sition (L=0.5); (C) pentacene on a 1 nm water on Si02; 
(D) pentacene/Si02 measured at the quarter-Bragg position 
(L=0.25). 



B. Layer-wise Rate Equation Model of Epitaxial 
Growth 



A(q z )<x [(po- Pl )e^ LT + (p 1 - P2 )] . (5) 



Thus, the scattering amplitude of a rough film, which 
corresponds to the center of the circular trajectories of 
ideal SF growth in the left-hand plots in Fig.[3j also corre- 
sponds to the scattering amplitude of the smooth, buried 
interfaces alone. 



Even if the scattering parameters described in section 
|III A| are precisely known, it is clear that the coverages 
9 n (t) at time t are not uniquely determined by a sin- 
gle intensity value I(t). To determine these coverages, 
additional information about the surface evolution is re- 
quired. The most common approach is to construct a 
deterministic, parameterized model for the evolution of 
coverages 9 n (t). Such models take the form of coupled, 
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differential equations: 



Fig. [2j, is written: 



d6n 

dt 



F(kj 



'n 1 "n+1 j • ■ 



(6) 



where the functions F represent the additional informa- 
tion imposed on the system, and may or may not depend 
explicitly on a subset of coverages and time t. The pa- 
rameters kj determine the 0„(i)'s, which are substituted 
into Eq. [2] to calculate I(t). 

A large variety of models of the form in Eq. [6] have 
been developed over the last several decades. Most of 
these originated and evolved in the context of specific 
experimental techniques, such as Auger Electron Spec- 
troscopy (AES) 33 36 , ion-beam assisted depositiorPBEH, 
and gJHEECp Bll Models used to analyze XRR 
dat d 13 * 14 * 16 * 23 * as described here have virtually all drawn 
directly from those of Ref. [40 . More recently, Trofi- 
mov et a?P2lED developed a new variant of this class of 
models, specifically for the purpose of linking atomistic 
kinetics to morphology in thin, multilayer films. 

It is convenient to categorize these models according 
to how much they attempt to incorporate atom-level ki- 
netics. Mod els by Kariotis et a/P^, Trofimov et a/EHMl^ 
and Koponer l 37 ! 38 ! draw explicitly on nucleation theorjEl 
and include additional equations to Eq. [6] representing 
the adatom and island densities on each layer. The next 
simplest models, most notably those introduced by Co- 
hen et al™, but including Refs. [35], [36], and [39], do 
not explicitly attempt to model atom-level kinetics, but 
approximate these kinetics in the form of Eq. [6] For 
example, the "distributed" model described in Ref. [ID] 
includes a mean-field representation of the step-density 
of a layer as a function of coverage 9, which in turn con- 
trols the amount of downhill interlayer diffusion. Finally, 
Braun et al^l have developed a version of Eq. § which 
attempts, explicitly, to avoid modeling atomistic kinetics. 
Rather, the equations Eq. [6] are coupled via intermediate 
functions J n (t), chosen to be as simple as possible con- 
sistent with producing physically reasonable results for 
the functions 9 n (t). 

Below, we describe implementations of three of the 
models described above: the empirical model described 
by Braun et al^, a simplified version of th e model intro- 
duced by Trofimov et aZp2l, and our variant! ^ ! 20 ! 48 ! of the 
distributed model introduced in Ref. [40 . In each case, 
we introduce analogous modifications that were found 
necessary to fit data obtained for our model systems, 
pentacene/SiC>2 and DIP/Si02- Specifically, successful 
fits, described in section |TV} require that the kinetic pa- 
rameters change with layer number for at least the first 
three layers. For the case of DlP/SiC^, it is additionally 
required that the sticking coefficient of incident molecu- 
lar species vary with coverage of the first layer. 

First, we present two modifications to the so-called 
distributed model of Ref. [40] that widen its range of 
applicability and connect, within the limits of the mean- 
field approach, to well-defined physical quantities. This 
modified Cohen (mC) model, (depicted schematically in 



dt 



l) 



where R n and a n are the net deposition rate from vapor 
and downward interlayer transport probability into the 
n layer. The parameters R n and a n implicitly represent 
all of the kinetic processes involving molecular attach- 
ment and transport, respectively. Rather than Cohen's 
original form for a„, we use 



) = 



k-d(6 n ) 



k-d{0 n ) + k + d{6 n+1 ) 

e- E 'd(6 n ) 
e- E 'd{6 n ) + d{6 n+i y 



(8) 



where d(9) represents the average step-edge density of a 
layer with coverage 9. The factors &-(+) represent the 
rate of molecular attachment at downhill (uphill) steps, 
and e = k-/k + . In the assumption of irreversible 
attachment, E' = E-^s/kT, where E-es is the Ehrlich- 
Schwoebel barrier for interlayer transport! 4 - . In general, 
the function d(9) could depend on layer number, for ex- 
ample if the island density changes as a function of layer 
thickness or island shape. Since such differences are vir- 
tually indistinguishable from differences in E' , we use 
the same function d(9) for each layer, but typically allow 
variations in E' as a function of n. 

The rate equations described by Eq.[6]do not explicitly 
prevent overhangs, i.e. solutions in which 9 n (t) > 9 n -\(t) 
for some n and t. Overhangs can occur, for example, if 
a n > «„_!, so that the interlayer transport into layer n is 
much greater than the transport of layer n into layer n— 1. 
In our implementation of the Cohen model, overhangs are 
prevented by forcing a n to approach as 9 n approaches 
9 n -i- 

The parameterization in Eq. [8] fixes a drawback of 
the original form, in which a„ is defined as^l a n = 
Ad(9 n ) I (d(9 n ) + d(9 n+ i)). In this case, A = 1 corre- 
sponds to E' = — 00 in our model, whereas < A < 1 
corresponds, approximately, to E' > 0. The region that 
is mathematically forbidden in Cohen's original form, 
— 00 < E' < 0, is precisely the region of parameter space 
where extended LBL oscillations are expected. As noted 
in Ref. [20J that model is unable to reproduce such be- 
havior. 

Our second modification to the Cohen model is to em- 
ploy recent work of Tomellini et concerning the form 
of d(9): the mean-field step-edge density of a growing 
film comprised of 2D islands. Two forms of d(9) are ob- 
tained, depending on whether or not islands rearrange 
when they merge, reducing their total perimeter. The 
zero and complete rearrangement limits are labeled "im- 
pingement" and "coalescence" regimes, respectively, and 
are found to have step densities d lm and d c ^^. 



d im (6) = 2^No{l - 



In 



(1-9) 



1/2 



(9) 



G 



d co {6) = ^9(l~6)cxp(- 



2(1- 



(10) 



Equation [TU] reflects the fact that if islands rearrange 
upon coalescence, the density of holes remaining in the 
surface near layer completion is smaller than the island 
density at nucleation. This implies that the step-edge 
density is also smaller, e.g. d(l — e) <C d(e). Real sys- 
tems are expected to exhibit behavior intermediate be- 
tween that of Eqs.[9|and|10[ In practice, we parameterize 
this form so that the degree of coalescence may be tuned 
by a parameter 5, < 5 < 1: 



d(9) = 5d co (9) + (l-S)d im (9). 



(11) 



In order to solve Eq. [7] using Eqs. |8]and[TT] it is neces- 
sary to "nucleate" layers 2 and above, since a n otherwise 
remains equal to one, precluding nucleation of layer n+1. 
This is done by computing d(9 n+ i) using a small value 
(e ~ 1CU 5 ) of 9 n+ i when 9 n exceeds a critical coverage 
9 n cl . Explicitly: 



1 



< 



a n (9 n ,9 n+1 + e) 



(12) 



Using Eqs. [8] and |11| results obtained by numerically 
solving Eq. [7] may be tuned, via the parameter £", from 
3D growth (E 1 > 1) to perfect LBL growth (£" < -1). 
More complicated behavior can result from allowing dif- 
ferent values of E' for different layers. For pentacene 
deposition on SiC>2, we find that reasonable results (dis- 
cussed further in section IV) are obtained by employing 
several values of E', e.g. E\, E 2 for the first two layers, 
followed by asymptotic approach to E 2 + AEn for layers 
n > 3, e.g. E n > 3 = E 2 + AE N x exp((n - 2)/N ). For 
DIP/SiC>2 and other systems exhibiting growth rate ac- 
celeration, Ri is allowed to differ from R n> i, simulating 
a difference in sticking coefficient for molecules incident 
on the growing film as opposed to the bare substrate. In 
summary, then, our model for real systems involves ten 
growth parameters, R\, R n >i, E\, E 2 , AE N , 9 l cr , N , 
9i, cr, #jv,cr, and 5. 

The second model in our comparison, introduced by 
Ref. [15] and referred to below as the Braun/Kaganer 
(BK) model, is explicitly simpler than those described in 
Ref. [10] • Equation [6] takes the form: 



d9n 
dt 



(13) 



where the intermediate functions J n depend explicitly on 
time as 



J ™=2 ,[ 



tanh 



t t n 

Jn 



(14) 



Here, the times t n are determined directly by the rates 
R n , while the parameters Bn , which are generally chosen 
to increase monotonicalljff^, have the effect of determin- 
ing the film roughness. As with our model above, we 



find that additional variation in the values of /3„ for the 
first few layers is required to obtain a good fit to data. 
We define (3i and f3 2 as free parameters, while for layers 
n > 3, f3 — fi 3 (n — 2) a . Because of the explicit time 
dependence of the intermediate functions J„, this model 
does not strictly allow for a coverage dependent growth 
rate. Nevertheless, allowing different values for R n does 
approximate such behavior. As with the mC model, the 
value of Ri is allowed to differ from that of R n >i for sys- 
tems exhibiting growth rate acceleration. To calculate 
the t n 's, we make the explicit assumption that layer N 
is completed at time t N = 1/R\ +X^=2 l/-^n>ij leading 
to t n = nR n> \ I \R\R n> \) . Thus, the BK model includes 
six parameters, R\, R n >i, Pi, j3 2l /?3, and a. 

The last model employed for our comparison is a sim- 
plified version of an atomistic, rate-equation model devel- 
oped by Trofimov et aZP^HUl. Unlike the layer-wise mod- 
els just described, this model includes three rate equa- 
tions for each layer, describing the rates of change of the 
adatom density, island density, and coverage. Downhill 
transport is controlled by a so-called "feeding zone" (FZ) 

of each layer n: atoms that are incident onto the feed- 
ing zone of layer n remain on that layer and hence con- 
tribute to layer n+1. Adatoms that land on top of layer 
n but outside this feeding zone diffuse downward, thus 
increasing 9 n . A clear advantage of this model is the abil- 
ity to directly connect morphology, including in-plane pa- 
rameters such as the nucleation density, evolves as a func- 
tion of layer number, and as a function of real, physical 
parameters such as D/J, the ratio of adatom diffusivity 
to the incident flux. For example, it is found that, even 
for homoepitaxy, the saturation island density decreases 
as a function of layer number during LBL growth^. This 
phenomena has recently been directly observed in SrTiC>3 
homoepitaxy using pulsed laser deposition^. 

For the limited problem of examining layer coverages, 
the FZ model may be re-parameterized so that it is rep- 
resented by just one equation per layer. Examination of 
Ref. [U] shows that all of the atomistic physics is con- 
tained in the evaluation of the critical coverage for each 
layer 9 nfix . Defining this as a fit parameter, we have: 



dK \R X {\-Qx)- 
dt 1 -R„>i(£,j_i 



R n >i(9i - £i) > n = 1 
■£n) , n > 1, 



(15) 



where 



'0 , 9 n < 



l_ e ~ [V- ln(l-0„)-V-ln(l-0n,cr)] ' 



> 



. (16) 

As in the models above, we have incorporated variation 
in the sticking coefficient of incident species through the 
parameters R\ and R n >i 7 and obtain reasonable fits by 
allowing 9 njCT to change with layer number. The analy- 
sis shown below incorporates six growth parameters, 
Rn>i, #i,crj 9 2 cl , #oo,cr7 and a parameter N determining 
the asymptotic approach to #oo cr . 



IV. COMPARISON TO EXPERIMENT 
A. Pentacene/Si02 

The AFM results in Fig. [TJ3-E obtained by hyperther- 
mal growth are consistent with those of prior . voluminous 
work on (thermal) pentacene/Si0 2 MBB 16 ' 29 ' 30 ™ !. 
Growth begins in a LBL mode, so that the first layer is 
nearly complete before the second layer nucleates. Subse- 
quently, roughness increases quickly, resulting in a late- 
time morphology characterized by large, multilayer is- 
lands. The precise nature of this transition varies with 
growth conditions — particularly temperature^!. 

Figure |4]A_ reproduces the data from the thickest film 
shown in Fig. [T] along with fits to the three models de- 
scribed in section [ill B| Figures [4J3-C show the evolution 
of layer coverages and rms roughness, respectively, re- 
sulting from the fits in Fig. [4]^, in addition to the rms 
roughness values (Fig. |4p) obtained from AFM on each 
of the four films. Error bars in Fig(4p, visible only for 
the thickest film, correspond to standard deviations of 
values obtained from different quadrants of the images in 
Fig. [TJ3-E. As discussed further below, however, as much 
as 15% variation in growth rate was subsequently found 
across the 1-cm wide substrate, all of which is sampled by 
the x-ray beam. Therefore, these error bars may under- 
estimate the variation in film morphology as observed by 
the x-ray measurements. Nevertheless, reasonably good 
agreement is also obtained in rms roughness evolution 
among the models and the AFM results, particularly at 
early times. 

The parameters for the fits in Figj4] are shown in Ta- 
ble HJ Since deposition was performed directly on clean 
SiC>2, as confirmed by x-ray reflectivity prior to growth, 
no interfacial layer was included (pi — 0, r = 0). For 
Po, a bulk mass density of 2.2 g/cm 3 was used, corre- 
sponding to the assumption that electron density is pro- 
portional to mass density. For this case, no growth rate 
acceleration is observed: allowing R n >i to vary from R\ 
does not improve the fit. Thus, R n >i = R\ for the fits 
shown. All models reproduce the conspicuous features 
of the data, giving R-squared values of 99.4% or higher. 
However, the fits are not perfect: the BKand FZ models 
give comparable quality fits, with xl values^ of 35 and 
34, respectively, indicating that the differences between 
the model and data are statistically significant!^. For 
the mC model, the best fit, with xt = 7, was obtained 
with Ai?2 and #oo,cr at 0, and the remaining parameters 
as shown in Table III Error bars in the table correspond 
to an increase in ^ by 1, after re-optimization of the 
remaining parameters. 

The fact that the xi values differ significantly from 1 
reflects, in part, the high accuracy of the data: the aver- 
age intensity I in Fig. |4|V is 2.4xl0 4 , so that the mean 

statistical uncertainty in the data (^s/l/IJ is 0.7%. Al- 
ternatively, we compare the mean absolute residual with 
the mean of the data (\I — I m \ /I m )- For the fits in Fig.|4j 
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FIG. 4: (Color online) a) (open circles) XRR data obtained 
near the pentacene anti-Bragg position during pentacene de- 
position on Si02- Also shown are best-fit curves to the mC 
(solid line), BK (dotted line), and FZ (dashed line) mod- 
els described in the text, b) Layer coverages from the three 
models shown in (a), c) Evolution of root-mean square rough- 
ness obtained from the three models shown in (a) and (b). 
Filled squares are measured rms roughness values obtained 
from the AFM images shown in Fig. [TJ3-E, corresponding to 
film thicknsesses of 0.34, 1.4, 2.5, and 4.4 ML. 



the mC model differs from the data by 1.4%, the BK and 
FZ models by 3.2%. 

In view of the greater number of parameters in the 
mC model, the better fit is not unexpected. The ques- 
tion remains whether this statistically improved fit cor- 
responds to more accurate layer coverages in Fig. [4j3. 
The rms roughnesses appear to suggest that this is the 
case, since the measured rms roughness at t = 500s cor- 
responds more closely to the mC model than to the oth- 
ers. A more detailed comparison of the models and films 
is given in Fig. [5j which shows height distributions ob- 
tained by the AFM data in Fig. [T] as well as the height 
distributions obtained from each of the best-fit simula- 
tions in Fig. [4] In each plot, discrete height distribu- 
tions (represented as black bars) were obtained by fitting 
the continuous distribution to a sum of several Gaussian 
peaks, one per layer, and equating the area of each peak 
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AE N 
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TABLE I: Fit parameters for the three fits shown in Fig. [4] 



to fractional, exposed occupancies for each layer c„. If 
the c„'s are normalized such that J2 n c « = 1> they are 
related to the layer coverages 9 n as: c„ = 6 n — 6 n +\. For 
AFM images, the indexing n of layers is chosen such that 
the total thickness = J2 n nc n — S« ®n most closely 
matches thickness estimates provided either by the x-ray 
data or, alternatively, by the sub-monolayer growth rate. 
Although not shown, we assume a maximum error of 15% 
in the c n 's, arising from the growth rate inhomogcncity 
across the film surface, as discussed above. Figure [5] re- 
veals that the agreement between the mC model (dark 
gray) and AFM is not conspicuously better, and may in 
fact be worse than that of the other simulations. Thus, 
the lower value of xt does not necessarily imply a more 
accurate representation of the true morphological evolu- 
tion of the film. 

The parameters in Table [I] warrant several comments. 
First, the values of p2 for all three fits are close to, but 
somewhat larger than that of bulk pentacene, 1.3 g/cm 3 . 
The sign of this discrepancy is consistent with the ex- 
pected error due to the uniform slab approximation. The 
fact that the electron density in pentacene is slightly con- 
centrated around the molecular center (along the c axis) 
has the effect, at L values below the first Bragg Peak, 
of increasing the scattering amplitude of this layer com- 
pared to a uniform slab with equivalent average density. 
A larger, uniform electron density compensates for this 
difference. Second, although the data are nominally ob- 
tained at L = 0.5 (the anti-Bragg position), good fits 
require L w 0.46 — 0.48, and an offset between the sub- 
strate and the first layer of A z — A x c « 3A. Possible 
contributions to the deviation of L from 0.5 include the 
experimental uncertainty, defined by the detector slits to 
be AL = 0.03, and a change, with layer thickness, of the 
pentacene d-spacing c. Such a change has been reported 
by Fritz et a/p2. The offset A z is assumed to correspond 
to a thin interfacial region at the SiC^/film interface. 
The remaining parameters in Table |TJ those describing 
the growth morphology, all show effectively the same, 
monotonic trend from smooth to rough growth. Specifi- 
cally, rougher growth corresponds to increasing values of 
/?„ in the BK model, decreasing values of 6* njCr in the FZ 
model, and a combination of increasing E n and decreas- 
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FIG. 5: Height distributions obtained from the AFM im- 
ages in Fig. [l] and from best-fit simulations shown in Fig.|4j\. 
Open circles represent height distributions obtained directly 
from the AFM data, solid lines indicate fits of these data 
with a sum of several Gaussian distributions (one per layer), 
and black bars indicate the area of each Gaussian distribution 
so obtained. Dark grey, light gray, and white bars represent 
height distributions obtained from the mC, BK, and FZ mod- 
els, respectively. 



ing #„, cr in the mC model. 

Finally, we note that the three models predict growth 
rates R x which vary from 0.0081 to 0.0089 ML/second. 
Evidently, even for the fairly straightforward problem of 
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determining the growth rate of this simple system, the de- 
tails of the model strongly affect the results. The origin 
of this variation is the link between the onset of rough- 
ness and the appearance of a peak. In general, extrema 
in the scattering data occur when the growth rate of one 
layer overtakes that of the layer below it. For perfect 
LBL growth, this moment coincides with layer comple- 
tion, producing cusp-like peaks. But for a film undergo- 
ing roughening, this transition can occur at lower cov- 
erage for each consecutive layer, resulting in a shorter 
time between peaks than the time to deposite 1 ML. If 
the roughening transition is abrupt, the growth rate can 
appear to accelerate. This is well illustrated by Fig. |3]\. 
Growth begins at t = 50s, and the first two local max- 
ima occur at t — 167s and t = 275s, corresponding to 
growth rate estimates of 0.0085 and 0.0093 ML/second. 
The three different models represented in Fig. |4]A_ model 
the onset of roughness in slightly different ways, resulting 
in three different estimates of growth rate. 

The variation in best-fit growth rates among different 
models, all of which fit the data well, raises the question 
of whether there is an independent means of measuring 
growth rate and, especially, of determining whether or 
not growth rate acceleration has occurred. We rei terate 
that for the hyperthermal growth method used her e! 19 * 25 *, 
the use of a quartz crystal monitor is not possible due 
both to the narrowly-directed beam profile and the pos- 
sibility that sticking coefficients are generally substrate- 
dependent. We present two alternatives. First, the cusp- 
like nature of the first local maximum in Fig. [4]A. sug- 
gests that this peak closely coincides with completion of 
the first layer. Thus, provided that we are confident that 
growth rate acceleration does not occur, and that the 
surface morphology does not develop large height asym- 
metries during the first layeJ^, the growth rate estimate 
given by the peak, 0.0085 ML/second should provide an 
accurate estimate of the growth rate. This growth rate 
falls within the results from the three fits. 

A second alternative measure of the growth rate can be 
made using the height distributions represented in Fig. [5] 
to obtain a plot of thickness vs. time. As noted above, if 
layers in an AFM image are correctly indexed, the thick- 
ness O is equivalent to the center-of-mass of the height 
distributions ^2 n nc n . The circles in Figure m\ represent 
the thicknesses of the four films in Fig. [T] obtained in 
this manner, along with the best-fit line to those data. 
Clearly, a line describes the data well. Unfortunately, 
the growth rate obtained in this fashion, 0.0096±0.0002 
ML/s, does not fall within the growth rates obtained from 
fits, and is not consistent with 0.0085 ML/s obtained by 
simple inspection of Fig. [4|A. It is difficult to reconcile the 
x-ray data with a growth rate of 0.0096 ML/s, since this 
growth rate would imply that the sharp, cusp-like peak 
at t = 167s occurs at a thickness greater than 1 ML. 

One possible source of error of the AFM-obtained 
thicknesses in Fig. [6] is the finite size of the AFM tip. 
The finite spatial resolution resulting from tip size will 
always result in an overestimate of the film thickness, 
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FIG. 6: Thickness, determined by AFM as described in the 
text, vs. time of (circles) pentacene films represented in Fig. 
[I] and (squares) DIP films represented in Fig. [7] Also shown 
are linear fits to the AFM-determined thicknesses. The best- 
fit line to the pentacene data intersects the origin, indicating 
an absence of growth rate acceleration, and a growth rate of 
0.0096 ± 0.0002 ML/s. For DIP, only the later three point 
are used for the fit, yielding a late-time growth rate of R„>i 
of 0.0156 ± 0.001 ML/s, and a negative intercept, indicating 
growth rate acceleration. 



since holes will appear smaller than their true size, while 
islands will appear larger than their true size. We esti- 
mated this error by finding the total perimeter of all of 
the islands for several layers and images in Fig.[l^-E, and 
multiplying by the estimated error, 4 nm, of the mea- 
sured, lateral position of the step due to a finite tip radius 
of 10 nm. This lateral position is assumed to correspond 
to half of the full 1.5nm height of a single ML step, and 
the AFM tip radius is obtained from SEM measurements 
of a nominally identical tip to that used to obtain the im- 
ages in Fig. [TJ3-E. The resulting error is approximately 
1% of the total image area, much less than the 10% dis- 
crepancy between the simulations and Fig. [6j Of course, 
because the islands are fractal, our measurement of is- 
land perimeter is an underestimate. A simple estimate 
of this error, considering the AFM tip size, molecule size, 
and estimated fractal dimension df ss 1.8p3 suggests our 
perimeter measurement to be in error by no more than a 
factor of 2. 

Another explanation for this discrepancy is that the 
growth rate is not uniform over the portion of the film 
measured by x-rays. Because these measurements are 
performed at an incidence angle of 1°, the 0.5mm x- 
ray beam probes the entire 10 mm width of the sub- 
strate. Subsequent AFM measurements (not shown) con- 
firm that the growth rate near the edge of the film is 15% 
smaller than that near the center. Thus, the average 
growth rate measured by the x-ray beam is smaller than 
that at the center of the film, where the measurements 
in Fig. [I] were obtained. 

Regardless of the disagreement in growth rate implied 
by the x-ray and AFM analysis, the two approaches agree 
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FIG. 7: AFM images of four thin films of DIP/Si0 2 of in- 
creasing thickness grown in immediate succession on four ar- 
eas of the same substrate: (A) 0.97 ML; (B) 2.06 ML; (C) 
4.17 ML; (D) 11.16 ML. All images are lO/ttmXlOjum 



with regard to the absence of growth rate acceleration. 
Referring again to Fig. [6j this conclusion comes from the 
fact that the best-fit line of the pentacene data intersects 
the origin to within 0.02 ML. This finding contrasts the 
case of DIP/Si02, also represented in Fig.[6]and discussed 
in detail below. 



B. DIP/Si0 2 

Figure [JJ shows AFM data obtained from a thickness 
series of DIP/Si02- As in Fig. [I] the films were grown 
in immediate succession on the same substrate, but at 
a substrate temperature of 89° C. This increased tem- 
perature may contribute to the two most conspicuous 
differences between Figs. [7] and [IJ namely, that DIP ex- 
hibits compact, rather than dendritic island morphology, 
and more persistent LBL growth. The later observation 
is clearly seen by comparing Fig. [4j\ to Fig. [8]A., which 
shows in situx-ray data obtained for the film shown in 
Fig.|7p. 

The squares in Fig. [6] represent thickness vs. time of 
the four films in Fig. u\ obtained by combining data in 
Fig. [JJ with qualitative analysis of in situx-r&y data ac- 
quired during growth. In contrast to the pentacene case, 
a best-fit line to these four points clearly has a negative 
intercept, indicating growth rate acceleration. Because 
of this, only the later three points are used for the fit in 
Fig. |6j yielding a late-time growth rate of 0.0156 ± 0.001 
ML/s. 



To determine of the degree of growth rate acceleration, 
defined as (i?„>i — i?i)/i?i, we generally rely on quan- 
titative analysis of XRR data. However, for the special 
case of LBL growth, this, like the growth rate, can be 
obtained directly from the data in Fig. [6j In this case, 
the rate of growth of the first layer is described by: 

d ^ = R l {i-e 1 ) + R n>l e 1 , (17) 

and the remaining layers grow according to 

^ = R n>1 (l-9 n ). (18) 

Solving this system yields the total thickness = @n 
as a function of time: 

f «i ( e (R n >i-Ri)t _ i\ t<u 

\l + Rn>i{t - h) , t>t u 

where h = \og((R n>1 - Ri)/Ri + l)/{R n >i - Ri). We 
now identify the linear regime in Eq. [19] with the best-fit 
line to the late-time thickness data in Fig. [6j Setting the 
intercept b of that line equal to 1 — R n> \ti yields 

(l-&)(/^-l)=log(/^), (20) 

which can be numerically solved for R\. For the DIP 
data in Fig. [6) R n>1 = 0.0156 ML/s and b = -0.417 ML, 
yielding the result i?i = 0.0074 ML/s. This corresponds 
to a growth rate acceleration of 111%. 

In addition to the x-ray data in Fig. [8]^, Fi g. [8] sh ows 
fits of these data to the three models in section [Til B[ the 
simulated layer coverages and rms roughnesses resulting 
from these fits, and rms roughnesses obtained from the 
AFM data in Fig. [7] In general, the modeled x-ray inten- 
sities in Fig. [HJ\ show reasonable, but noticeably worse 
agreement than in Fig. [4]A_, and yield correspondingly 
worse xt values of 445, 756, and 509 for the mC, BK, 
and FZ models, respectively. As with the fits in Fig. [4j\, 
we also compare the mean residual to the mean value 
of the data (\I — I m \ //), obtaining values of 9.2%, 12%, 
and 10%. The reasons for the difference in fit quality 
between Figs.[4j\ and[8]A. are unclear. However, previous 
work on DIP /SiOj^ has demonstrated a clear structural 
change in a DIP thin film during growth of the first five 
layers. Such a change would have the effect of making 
some of the x-ray parameters discussed in section |III A| 
time-dependent. 

Referring to Fig. [8p, all three models compare very 
well with the AFM data up to t = 400s, corresponding to 
a thickness of 4 ML, beyond which the rms roughnesses 
predicted by the models diverge both from each other 
and from the actual film. At t = 800s, corresponding 
to a thickness of 11.2 ML, the mC and BK overestimate 
the roughness, whereas the FZ model underestimates it. 
These observations are reflected in more detail in Fig. |9j 
which compares the actual height distributions obtained 
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FIG. 8: (Color online) As in Fig. [I], a): (open circles) XRR 
anti-Bragg data during DIP deposition on Si02, along with 
best-fit curves from the mC (solid line) , BK (dotted line), 
and FZ (dashed line) models, b) Layer coverages from the 
three models shown in (a) . c) Evolution of root-mean square 
roughness obtained from the three models shown in (a) and 
(b), with values obtained from ex situ AFM shown as filled 
squares. 



both from the AFM data and from the models. The cor- 
respondence between the simulated and measured height 
distributions for the three thinnest films is very good, 
whereas all three simulations depart significantly from 
the measured height distribution of the thickest film in 
the series. 

Echoing the analysis of pentacene/Si02, the mC 
model, which yields the lowest value of xt compared to 
the other two, nevertheless does not provide a conspic- 
uously better representation of film growth. Moreover, 
both the mC and BK models tended towards second, 
deeper xt minima (xt =285, 417, respectively) charac- 
terized by the disappearance of the shallow local maxi- 
mum near t — 150s (see Fig. 10 <Y). As demonstrated by 
Fig. [Top, the film morphology implied by these alterna- 
tive fits are inconsistent with the AFM data. To avoid 
these minima, only a subset of parameters could be al- 
lowed to vary simultaneously, taking care that the best-fit 
model reproduced, at least weakly, all of the extrema in 
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FIG. 9: Height distributions obtained from the AFM im- 
ages in Fig.[7]and from best-fit simulations shown in Fig.[8]A 
Open circles represent height distributions obtained directly 
from the AFM data, solid lines indicate fits of these data 
with a sum of several Gaussian distributions (one per layer), 
and black bars indicate the area of each Gaussian distribution 
so obtained. Dark grey, light gray, and white bars represent 
height distributions obtained from the mC, BK, and FZ mod- 
els, respectively. 



the data. 

The parameters for the fits in Fig. are shown in 
Table |Tl} As for the analysis of Fig. [4j\, these fits presup- 
posed the absence of an interfacial layer (pi =0, r — 0), 
and employed po = 2.2g/cm 3 . For the mC model, #i. cr , 
02.cr, and #oo jC r were all fixed at 0. 

An interesting feature of Fig. |8p, in both the AFM 
data and the models, is that the rms roughness at film 
thickness O — 2 is smaller than that at = 1. This is 
reflected in the fact (see Table [TTJ) that the parameters 
controlling interlayer transport for all three models do 
not change monotonically, as in Table [TJ The cause and 
meaning of this behavior are outside the scope of this 
paper: however we suggest that it could be related to the 
structural evolution of the DIP filrrPH r the dependence 
of growth rate on coverage. 

The key difference between Tables [n] and [TJ is the ad- 
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TABLE II: Fit parameters for the three fits shown in Fig. [5] 



ditional parameter, R n >i, corresponding to the late-time 
growth rate of the film. Above, we found that the AFM 
data in Fig. [6] gave a late-time growth rate of 0.0156 
ML/s and a growth rate acceleration of 111%. Refer- 
ring to Tab. |TlJ we find that this late-time growth rate 
agrees well with that obtained from all three models. The 
growth rate acceleration for the mC and FZ models yields 
112% and 122%, respectively, in excellent agreement with 
the result from AFM. The BK model, which we reiterate 
is not intended to model growth rate acceleration accu- 
rately, gives a more modest growth rate acceleration of 
56%. We conclude that, at least for systems exhibiting 
relatively smooth growth at early times, the mC and FZ 
models provide an accurate means of determining growth 
rate acceleration. 



V. SUMMARY AND CONCLUSIONS 

In this paper, we have provided the first comparative 
study of the problem of obtaining surface morphology 
from in situ, time resolved XRR data. We described 
three mean-field models of thin film growth, along with 
detailed examples of their application to quantitative 
analysis of XRR data obtained at the anti-Bragg posi- 
tion during growth. Two sets of organic thin films were 
grown using hyperthermal deposition, pentacene/Si02 
and DIP/Si02- For each system, in situXKR data was 
obtained during growth of four films, grown to different 
thicknesses, under nominally identical conditions. The 
XRR data were fit to each of the the three models, result- 
ing in detailed simulations of the time-dependent mor- 
phology of each film. Finally, these simulated morpholo- 
gies were directly compared with AFM data from each 
of the four films to critically evaluate the quality of each 
simulation. 

We find, first, that all of the models provide good de- 
scriptions of both the XRR data and, at early times, the 
surface morphology of an evolving film. This is signifi- 
cant, since the evolution of layer coverages incorporates 
the rate of interlayer transport, which controls the onset 
of roughness. Second, we find that the model which fits 
a particular data set best statistically, i.e. which gives 



the lowest value of xti does not necessarily provide the 
most faithful reproduction of the true surface morphol- 
ogy. The models we used all incorporate analogous pa- 
rameters determining interlayer transport as a function 
of layer coverage, but these parameters represent fine 
distinctions regarding intralayer atomistic kinetics and 
morphology. The fact that better fits do not correlate to 
improved agreement with the true surface morphology 
suggest that our measurements are not sensitive to such 
fine distinctions. 

We understand this lack of sensitivity to fine detail as 
resulting from the limited amount of information con- 
tained in the intensity at a single point q z on the reflec- 
tivity curve. As stated above, in situstudies performed 
at multiple q z values have been demonstrated, both in 
angular dispersive^ and energy dispersive modes using 
white bearrp2. However, these techniques both suffer 
from poorer time resolution than in single point mode. 
In angular dispersive mode, this results from the require- 
ment of moving the sample and/or sample chamber dur- 
ing the growth. In energy-dispersive mode, the time res- 
olution is limited by detector dynamic rang tP 1 , and the 
fact that the signal is so much stronger near Bragg peaks 
compared to the anti-Bragg position. A variation to this 
approach would be to combine an energy-dispersive de- 
tector with a finite bandwidth beam, for example 10%, 



obtained with wide bandwidth multilayers^ in conjunc- 
tion with a bend magnet or wiggler. Using a 10±0.5 
keV beam, an energy resolution of 0.2 keV would result 
in approximately 4-5 independent points on the reflec- 
tivity curve being obtained simultaneously, resulting in 
significantly more information with minimal loss in time 
resolution. 

Concerning the determination of growth rate: in the 
case of pentacene/Si02, which exhibits an abrupt rough- 
ening transition near 2 ML, different models obtain val- 
ues for the total growth rate that vary by as much as 
10%. For DIP/Si0 2 , which exhibit extended LBL oscil- 
lations, we find that fits of XRR data to the mC and FZ 
models both correctly extract the degree of growth rate 
acceleration. 

A key feature of the mC and FZ models is that they 
both incorporate measurable quantities, such as the layer 
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FIG. 10: (Color online) As in Fig. [4] , a): (open circles) The 
same XRR data as in Fig. [8] shown with alternative fits to 
the mC (solid line) and BK (dotted line) models, b) Layer 
coverages from the two models shown in (a), c) Evolution of 
rms roughness obtained from the two models shown in (a) and 
(b), with values obtained from ex situ AFM (filled squares 
These fits have lower xt values than those shown in Fig 
but rms roughnesses which are in clear disagreement with the 
AFM results. 



Appendix A: Thin Film Scattering Amplitudes 
Parameterized for In situXRR 



Here, we derive Eq. [2] from Eq. [I] resulting in approxi- 
mate values of Afii m /A su b and (/> su b for the technologically 
relevant case of an organic thin film on a substrate with 
a thin interfacial layer. Referring to Fig. [2j we begin 
by dividing the integral in Eq. [T] into separate volumes, 
treating the buried substrate, the interfacial layer, and 
each crystalline layer of the film as distinct. Since we are 
only interested in the specular intensity, we ignore the 
fact that the film layers are only partially complete, and 
instead treat the coverage 9 n as a modification to the 
density p 2 . This step ignores the qn contribution of the 
diffuse scattering } 51 * 57 !, which depends on both details of 
the system as well as instrumental resolution. Neglect- 
ing the static roughness of each interface, (which can be 
easily incorporated^ if desired) , and allowing only 9\ to 
be nonzero, Eq. [T] becomes 



A(q z )/A = p 



/ dze~ iq ' z + pi dze 

J-oo J-Ti 



iq z z 



c+2A 2 



dzp 2 {z)e- iq >\ (Al) 



where A is the sample area, and 2A Z is a correction term 
for the height of the first layer. 

We next carry out the integrals, also allowing layers 
ri > 1 to be nonzero: 



A(q z )/A = -[p 1 + (p a -p i y^)] 

1z 

N 

+ A mol (q z )e-^ c / 2+A ^ ]T 9 n e-« n -V«* c . (A2) 

ra=l 

The term A mo i(q z ) is the molecular structure factor, de- 
fined as 



A mol (q z ) =J2e- iq * U J dz Pm (z)e-^ z , 



(A3) 



step-edge density and the critical coverage for next-layer 
nucleation, as fit parameters. Future work will aim 
to test whether these parameters, rather than only the 
layer coverages, can be extracted from XRR data during 
growth. For example, we will investigate how accurately 
the critical coverage parameter, # cr extracted from fits 
of XRR data to the FZ model, corresponds to indepen- 
dent measurements of 6 cr . We believe that such work 
will contribute significantly to the technologically impor- 
tant challenge of understanding, and ultimately control- 
ling, physical and experimental factors determining sur- 
face morphology in thin films. 



where the sum is over atoms in a unit cell, and £ m is the 
height of each atom relative to the center of the molecule 
(defined as z = c/2 + A z ) above the substrate, and p m {z) 
is the linearly-projected electron density of ato m m . Fi- 
nally, we factor out the term e _l9 ~( c / 2+Az ' in Eq. 



A2 



and 

convert to reciprocal lattice units L = q z /(2Tr/c), defin- 
ing t = Xi/c, A = A z /c. Equation A2 then becomes 



•11 '7-. )/ -X - ^T/e <7riC1+2A) (/>i + (Pa - Piy 27rL l 



N 



A moX {q z )Y,0ne- Kn - 1)2 *\ (A4) 



in which the first term is equivalent to A SVL be^ uh in Eq 
Equation | A4| thus accomplishes the goal of explicitly 
taining the parameters 4> su b and Afii, n /^4 su b in Eq. [2] in 
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terms of physical parameters of the film. We note too, 
that for centrosymmetric molecules such as pentacene, 
is pure real, so that at the anti-Bragg position 
(L = 0.5), the second term contributes only pure real 
terms to the complex sum. For non-ccntrosymmetric 
molecules, A z may still be chosen such that A mo \(q z ) is 
pure real, in which case z = c/2 + A z may not correspond 
to the molecular center. Finally, we note that Eq. |A4| 



with explicit computation of A mo \(q z ) from Eq A3 can be 



used to calculate the entire specular reflectivity, includ- 
ing Bragg peaks, for reasonably thin films, e.g. T mm <C 
1/j.m. 

For large, organic molecules, the anti-Bragg position 
probes density fluctuations on an approximate length 
scale I = 2n/q z = c/L, which is large compared to in- 
teratomic distances. As a result, a molecular layer is 
well approximated by a uniform density slab. In this 
approximation, A mo \ can be written 



*mol 



2p2 2^L sin ( 7FjL )' 



(A5) 



We note that the simple case of homoepitaxy can be re- 
covered by setting the densities equal and t = 0. The 
uniform slab approximation at the anti-Bragg position 
L = 0.5 then gives the familiar resullpfl 



Jab oc |A A b| 



n 



i-2j]9 B (-ir 



(A6) 



Appendix B: Scattering Amplitude Evolution 
During LBL and SF Growth on Vicinal Surfaces. 

To compare trajectories of the total scattering ampli- 
tude during LBL and SF growth, we imagine an ideal, 
vicinal substrate with step-height c and terrace width 
W , and with the positive x direction perpendicular to, 
and pointing towards uphill steps. We choose the origin 
to coincide with the center of a terrace in x and height z 
of the first growing layer. The scattering amplitude A vic 
of a film with a total of Nl partially complete layers 
(compare with Eq. [2]) is then 



N W 



-i{q z c+q x W)k , 



^vic — ^ ' 
k=-Nv, 

A auh e^ + ]T e-^< n -^A ter (9 n ) 



, (Bl) 



where At eT ^ n (9 n ) is the scattering amplitude of layer n 
with fractional coverage 9 n . The left-hand sum in Eq. |Bl| 
is over Nt = 2Nw + 1 terraces, and defines the specular 
condition as q z c = —q x W. In perfect SF growth, each 
layer n commences only when layer n — 1 is complete, 
and grows from the up-hill step at x = W/2, advancing 
from right to left until reaching the down-hill step at 



-W/2 at 6 n = 1. Thus, 



W/2 



AterMOn) = f / die"**-*, (B2) 

Jw/2-e n w 

where / is the electron density per unit length of the 
growing layer. Along the specular rod, we can substi tute 



Bl 



— q z c with q x W in all of the terms e~ vlxC ^ n ~ 1 ^ in Eq 
Next we can bring these prefactors into the integrands 
and substitute variables, so that the right-hand term in 
brackets becomes 



E< 



-iq z c(n— 1) 



A 



-tor l«« i 



N L -2 „W/2~nW 

o = / E / dxe ~ iq " x 

n=0 J-W/2-nW 
pW/2-{N L -l)W 

+ f dxe~ lq * x (B3) 

J W/2-(N L -1+Gn l )W 

All of the integrals in Eq. |B3| may now be combined. 
Identifying the film thickness 6 = Nl — 1 + 9n l , we have 



n=l 



~iq z c(n—l) 



Ay 



9n) = f 



W/2 



dxe 



-iq :l ,x 



w/2-we 



J_ e -iq x W/2 

iq x 



[ e iq*We _ jj _ ( B4 ) 



for the total film scattering amplitude along the specular 
rod. Eq. B4 shows that the scattering amplitude in SF 
growth completes a circle in the complex plane whenever 
q x WQ = 27m, or equivalently (recalling —q x W — q z c = 
2irL), whenever = n/L. 

In LBL growth, deposited material lands randomly on 
each terrace, so the scattering amplitude of a partially 
filled layer is just Eq. |B2| evaluated at 9 = 1, then multi- 
plied by the actual coverage 9 n : 

r sm(q x W/2) 



A 



tor, LBL \ 



= 9 n fW- 



(B5) 



q x W/2 • 

Just like Eq. [2j Eq. |BT] with the substitution of Eq. [B5] 
results in regular, polygonal trajectories in the total scat- 
tering amplitude. 

The effect of considering a vicinal rather than singu- 
lar substrate is to reduce the scattering amplitude in the 
specular direction by a factor sm(q x W/2)/(q x W/2). At 
the anti-Bragg position, q x W — n for a vicinal surface, 
but q x = for a singular surface, implying that the scat- 
tering amplitude of a vicinal terrace is 2/ir that of the 
equivalent area on a singular surface. This decrease in 
intensity from the specular, anti-Bragg position is redis- 
tributed into off-specular angles due to diffraction from 
the steps, which form a blazed grating. Scattered inten- 
sity will appear at q x W = —it + 2im for all integers n. 
The total scattered intensity in the plane q z = 7r/c can 
be summed by adding up the square of Eq. |B5| over all 
allowed q x . Noting the identity (Eq. 1.422.4, p. 44), 



E( %\Tl{x + Tin) \ , 
— I = sur x 

V x + im 

71 — — CO V 



1 



OC 

y 

^— ' (x + 7m V 

n= — oo 



= 1, 



(B6) 
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we find that the total scattered intensity in the plane q z = 
7r/c for both vicinal and singular surfaces is {9 n fW) 2 . 
This is actually an example of a more general sum rule 
easily derived from Ref. the total scattered inten- 

sity at a given q z is independent of the surface morphol- 
ogy, provided that the lateral length scale of roughness 
is small compared to the coherence length. 

In homoepitaxy, the circular trajectory in Eq. |B4| for 
the scattering amplitude in the specular condition must 
result in c onstant intensity. In that case, and making use 



where A tcr is Eq. B5 evaluated at 6 
again used q z c 



-q x W. Adding Eq. B7 to Eq. B4 



1, and we have 
trans- 



of Eq. 



B5 



the term A sub e^ sub is evaluated as 



= 2/ 



1 — e iqiC 
sin{q x W/2) 



a iq z c/2 



(B7) 



lates the circle to the origin, so that the scattered inten- 
sity is constant with increasing 0. 
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